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We study the magnetic excitation spectra of Cu L2,3-edge resonant inelastic x-ray scattering 
(RIXS) from the spin-liquid ground state in two-leg spin ladder cuprates. Applying the projection 
method developed by the present authors, we have derived the formulas of the magnetic RIXS 
spectra, which are expressed by one- and two-spin correlation functions in the polarization changing 
and preserving channels, respectively. The one-spin correlation function includes both one- and 
two-triplon excitations and they are picked up separately by choosing rung wave vector — q a = n 
and 0, respectively. An application to Sri4Cu2404i reveals that the calculated RIXS spectrum 
captures well the dispersive behavior shown by the lower boundary of two-triplon continuum. By 
adjusting the geometrical configuration of the measurement, one-triplon dispersion around the zone 
center could be detectable in the RIXS experiment. The observed weak intensity in the higher 
energy region might be attributed to the two-spin correlation function, which could be detected 
more clearly at the M3-edge RIXS spectra. 



PACS numbers: 78.70.Ck, 72.10.Di, 78.20.Bh, 74.72. Cj 



I. INTRODUCTION 

Resonant inelastic x-ray scattering (RIXS) provides us 
with one of the rare opportunities to investigate various 
excitations in solids including charge, orbital, and spin 
degrees of freedoms:- Among them, extracting the mag- 
netic excitation by RIXS is rather challenging because 
to distinguish its intensity from other contributions such 
as the elastic signal and phonon excitations is hard to 
achieve due to the small energy transfer below 1 eV as- 
sociated with the magnetic excitation. However, it is 
worthwhile to pursue potential of RIXS as a promising 
probe for magnetic excitation, since it can survey a wide 
range of momentum transfer in the Brillouin zone and 
requires only small sample volume. These features con- 
vince us RIXS to be a choice complementary to the con- 
ventional inelastic neutron scattering (INS). 

Recently, progress in the experimental technique of 
RIXS enables us to carry out higher energy resolution 
measurement)^ for instance, in several transition metal 
oxides like the cuprates4^— Hill et al. have succeeded 
in detecting the magnetic excitation spectra peaked 
around 500 meV energy region at the Cu -RT-edge in 
La2Cu04 4. The following theories have revealed that the 
obtained signals are attributed to the two-magnon exci- 
tation from the antiferromagnetic (AFM) ground state 
brought about by the presence of the core hole potential 
during the intermediate states^— Since then, plenty of 
measurements have been reported out at the Cu i\-edge 
as well as Cu L 3 -edgeH-~— 

When we turn our attention to the L-edge RIXS, the 
situation becomes a little complicated. In the Cu L- 
edge RIXS process, the transition is between the 2p core 
and 3d states. The photo-excited electron eliminates the 



3d hole leading the 3d state to the closed shell. The 
absence of the spin degree of freedom at the core-hole 
site is similar to the situation of the non-magnetic im- 
purity problem in the spin system.! 16 ' 17 Constructing a 
relevant theory to handle such difficult situation is chal- 
lenging but attracting. Within the fast collision approx- 
imation, the momentum and polarization dependences 
of the magnetic excitation spectra at the Cu Z^^-edges 
were investigated^^— In another attempts, resonant en- 
ergy dependence of the RIXS spectra is described by the 
scattering operator inferred from elastic scattering^ 21 ' 22 

In order to extract the magnetic excitation from the 
closed 3d state at the core hole site in the intermedi- 
ate state, we have developed an effective theory to in- 
vestigate the magnetic RIXS spectrum at the Cu L 2 ,z 
edges i 23 ' 24 The theory treats the polarization and energy 
dependences of the spectrum faithfully, projecting the 
final states onto the possible spin excited states, from 
which the RIXS spectra are expressed as the form of the 
spin correlation functions. They consist of the two-spin 
as well as one-spin correlation functions. The latter is 
also found in the theory of INS, but the former is specific 
to RIXS. Even for the one-spin correlation function, its 
transferred momentum dependence turns out to be com- 
pletely different from the one in the INS theory reflecting 
the inclusion of the non-local magnetic excitation. 

We have applied our theory to two-dimensional system 
having the AFM ground state and confirmed the quan- 
titative effectiveness of the theory^ It has reproduced 
well the experimental data reported by Guarise et al. in 
S^CuO^C^. 13 Contributions of the multi-magnon exci- 
tations have been evaluated too. Then, we have applied 
our theory to one-dimensional Heisenberg chain having 
the spin-singlet ground state 24 that preserves rotational 
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FIG. 1: (Color online) Schematic sketch of the material and 
model, (a) Cu-0 network of ladder in Sri4Cu2404i • Filled 
(blue) and open (green) circles denote Cu and O sites, respec- 
tively, (b) Model of two-leg ladder 2x8 spins used to evaluate 



/m (wi)'s and fl (wi)'s. The periodic boundary conditions 
are adopted. The spin at site on leg 1 (encircled by dotted 
line) is annihilated in the intermediate state. 



invariance in the spin space. The RIXS spectrum, which 
can cover the magnetic excitations generated not only 
at the core-hole site but also around the core-hole site, 
has been derived in a way that manifests the rotational 
invariance. The results have been consistent with those 
obtained by other theoretical approaches £2r— Our result 
also suggests that contribution from the two-spin correla- 
tion function is detectable at the higher energy region in 
the er-polarization measurement. Unfortunately, at that 
time, no experimental data was available to compare with 
our results. 

In this work, we extend our theory to exploit the Cu 
^2,3-edge magnetic excitation spectrum in RIXS for two- 
leg spin ladder cuprate where the ground state is the 
spin-liquid retaining the spin rotational invariance. The 
low energy sector of the spin excitations are known as 
one- and two-triplon excitations^ Our analysis reveals 
that the RIXS spectra can trace the one- and two-triplon 
dispersions as a function of the transferred momentum 
through the one-spin correlation function. 

Then, we apply our theory to interpret the experimen- 
tal results reported by Schlappa et al. in Sri4Cu2404i.— 
This material is considered an experimental realization 
of a two-leg ladder structure (Fig. Q] (a))j2£ and has 
attracted of particular interest since Ca doped systems 



Sri4_ a; Ca a ;Cu2404i have exhibited superconducting state 
under pressure for x ranging from 11.5 to 13.6.— In 
the RIXS measurement at the Cu L3-edge (930.6 eV), 
Schlappa et al. observed the dispersive signals attributed 
to the lower boundary of two-triplon excitation contin- 
uum. Our theory succeeds in reproducing the momentum 
dependence of the RIXS spectral profile. In addition, we 
find one-triplon dispersion can be detected by rotating 
the sample around the 6-axis. The observed weak inten- 
sity in the higher energy region might be the contribution 
from the two-spin correlation function, mainly originated 
from two-triplon excitations with total spin S = 0. We 
argue that the contribution from the pure two-spin cor- 
relation function is distinguishable more clearly at the 
M3-edge measurement. 

The present paper is organized as follows. In Sec. [TTl 
we extend the theory developed in our previous papers 
aiming at the application to the Cu £2,3 edges in two-leg 
ladder systems. The RIXS spectra are expressed in terms 
of spin-correlation functions. In Sec. IIII1 the amplitudes 
leading to the spin-correlation functions are evaluated 
on a finite-size two-leg ladder cluster. An application to 
Sri4Cu2404i is shown in Sec. |lVl Section [V] is devoted 
to summary and discussion. 



II. THEORETICAL FRAMEWORK 

A. Initial state and magnetic excitations 

For the purpose of application to the two-leg spin lad- 
der system Sri4Cu2404i, we consider the system is at 
half-filling and is undopedM- with each Cu atom having 



one hole per site in the x 



orbital. 32,33 The x and 



y axes are defined along the Cu-0 bonds parallel to the 
crystallographic c and a axes, respectively, while the z 
along the b axis. Here, the a and c axes are along the 
rung and leg directions of the ladders. The ground state 
and the low-energy spin excitations are described by the 
5=1/2 antifcrromagnetic Heisenberg Hamiltonian with 
an additional four-spin exchange terms^ 



Hrr 



J\\ J]] J]] Si :T • S,-+i )T + J_l Si jT • Si,T 



i t = 1,2 i 

~\~ Jcyc 5 ' (^i,T ' S^_f.x,T) (^i,r ' Sj-j-i^) 
plaqucttc 

+ Jcyc ^ (Si,r ' S,;.y)(S^i iT • Sj+l^) 
plaqucttc 

~~ "^cyc ^ ' (^i,T ' Sj-)-l 1 r)(Sj-|-l ) r ' Siyf), (2.1) 
plaqucttc 

where the index i refers to the rungs. The symbol r 
G {1,2} discriminates legs, and r denotes 2 for t = 1, 
and vice versa. The exchange coupling constants along 
the legs and rungs are denoted as J11 and J±, respec- 
tively (Fig. [T](b)). In addition, four-spin coupling J cyc 
is included. 
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We briefly summarize the ground state and low energy 
sector of the spin excitations of H mag . The ground state 
is known to be the gapped spin-liquid due to the quantum 
fluctuation^ - — Unlike the AFM ordered state, the spin- 
liquid state preserves the rotational invariance. First, let 
us assume J|| = J cyc = 0. Since each rung is independent, 
the ground state is constructed by N/2 pairs of rung 
singlet where N denotes the number of spins. This state 
is called as rung singlet. The lowest energy excited state 
is realized by changing one of the singlets into triplet with 
the excitation energy, or equivalently a spin gap, J±. As 
the number of triplets increases, the excited states are 
called as one-, two-, three-triplon excitations, and so on. 

When J|| is turned on, this kind of simple picture may 
be invalid. However, it is known that even when Jm — > oo, 
the ground state is adiabatically connected to the rung 
singlet. 4 - Thus, classifying the magnetic excitation by the 
number of triplons gives a good description and the spin 
gap remains finitei 28 i 35 ~ — The presence of J cyc does not 
alter the nature of the ground and low energy excited 
states qualitatively for J cyc /J± ~ 0.25, although compli- 
cation sets in when the ratio grows, which is beyond the 
present interest^ 



B. RIXS spectra 

We extend our RIXS theory developed for spin- 
singlet one-dimensional Heisenberg chain 2 ^ to quasi one- 
dimensional two-leg spin ladder system showing the spin- 
liquid ground state. In the following, a concise version 
of the explanation of the theory is displayed, relegating 
a detail to ref. [24]. 

In the electric dipole (-B1) transition at the transition- 
metal L2 i 3-edge, the transition progresses between the 
2p-core state and the 3d state. This process, described 
by the electron-photon interaction Hamiltonian H mil is 
mediated by absorbing the incident photon with wave 
vector q i and energy Wj, and then, emitting the scattered 
photon with wave vector and energy u)f. The RIXS 
spectra may be expressed by the second-order El allowed 
process; J i 



stands for the lifetime broadening width of the core-hole 
and e n is the energy of the spin part in the intermediate 
state. The e n will be evaluated by the Hamiltonian -ff„ lag , 
constructed from -ff mag by eliminating the spin degree of 
freedom at the central core-hole site. 

In the intermediate state, the spin degree of freedom is 
lost at the core-hole site. The final state experienced such 
intermediate state may be expressed by So,i|g), So,2 1 g) , 
Si, ••• in the polarization changing channel, and 
So.i • Si t i|g), So,i • So,2|g), • • • in the polarization pre- 
serving channel, where suffix (0, 1) indicates the core- 
hole site. Paying attention to the non-orthogonality of 
these state, we project the final state on these state in 
the same way as carried out in our previous study on the 
one-dimensional system. We obtain 



E 



H int \n)(n\H int \^) 



core c n 



iT 



+/ 2 (1) (o; i )(Si,i + S_i,i) + /3 W (a; i )So,2 
+/i 1) (o; i )(Si, 2 + S_i, 2 )]|g}, 

^a f± ■ a l± [fi 2 \u t ) + /i 2) (^)S ,i • (S M + S_ M ) 
+ / 3 (2) (^)So,i • S , 2 + /i 2) (wj)S ,i ■ (Si, 2 + S_ 1)2 )] 

(2.3) 



where and a. fj_ stand for the polarization vectors 
of the incident and scattered photons, respectively, pro- 
jected onto the xy plane. The fji (cjj)'s are the coeffi- 
cients to be determined. They could be accurately eval- 
uated in a system having rather small size, since the rel- 
evant excited states are restricted around the core-hole 
site. 

Accordingly, the RIXS spectra for the polarizations 
a i{f) = ("?(/)> >"*(/)) are ex P ressed as 



W(q f a f ;q l a t ) = 2?r^ 



E 



(^ r \H int \n)(n\H int \^i) 



E a 



U!i — E n 



x d(E g + cj s ; - Efi - Wf), 



(2.2) 



with qi = (q l ,uj i ), q f = (q/,w/), |$;} = c\ i<Xi \g), |$/<) = 
Cq f a f \f), where \g) and \ f) represent the ground state 
and excited states of the matter with energy E g and Eft, 
respectively. The polarization directions of the incident 
and scattered photons are on and a/, respectively. The 
annihilation (creation) operator of photon with momen- 
tum q and polarization a is denoted as c qa (cl a ). The 
intermediate state is denoted as \n) with energy E n in 
the presence of the core-hole. Since E n includes the core- 
hole energy e core , we express E n = e coxe — iT+e n where T 



W{qftxf,q i OL i ) oc 
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af + Offx? 
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Y^ l \u3i;q c ,q a ,u) 



(2.4) 



where the first and second terms represent the con- 
tributions from the polarization changing and preserv- 
ing channels, respectively. The Y^(u)i;q c ,q a ,uj) and 
Y^> (wi\ q c , q a , w) are Fourier transforms of the one-spin 
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and two-spin correlation functions defined by 

F (1) (tJi;g c ,g Q ,cj) 
= J (z( 1 ) t ( Wj ;g c ,g Q ,0^ (1) (^;g c ,ga,0))e^dt(2.5) 

= J (Z( 2 )t( Wi ; 9c ,g Q ,t)z( 2 )(^; 9c ,g Q ,0))e^di(2.6) 

with 

Z W (oJi;q c ,q a ) 
_ e — <gerj-tg n (r-l)o r 



x [/i^( Wi )5^ T + /i lj ( Wi )(S| +1)T + s;_ 1>T ) 
+ / 3 (1) (a; i )5; iT + /i 1) ( Wi )(5; +liT + S3 ! _ liT )], (2.7) 



Z {2) (ui;q c ,qa) 



-iq c Tj -\-q a (r— l)a r 



x [/ 2 (2) (wi)(S i+1 ,T + S,_!, r ) + /| 2) (^)S i>7 



+ /4 ''( w i)(S :! -+l,T + Sj-iy 



S 7 - r . 



(2.8) 



The coordinate of the site at the j-th rung and r-th leg 
is denoted as (rj, (r — l)a r ) in the ca plane. The symbols 
q c and q a represent the wave numbers along the leg and 
rung directions, respectively. They are the corresponding 
components of vector , which is the projection of the 
scattering vector q = onto the ca-plane. In the 

two-leg ladder configuration, the latter takes two relevant 
values and ir/a r . Hereafter, the momenta q c and q a are 
measured in units of 1/cl and l/a r , respectively, when 
their numerical values are mentioned. 

Note that in the above expressions (|2.7[) and (|2.8I) . 
the magnetic excitations included are those at the neigh- 
boring sites within the two adjacent plaquettes linked 
to the central site as well as at the central core-hole 
site. Notice also that if the fast collision approxima- 
tion is adopted, Z^'fai] q c , q a ) and Z^ 2 ' (cJj; q c , q a ) be- 
come /i 1) (w ! )E J , T e- !,cr '-* (T - 1) ' lr S| T and 0, respec- 
tively, meaning that only the excitation at the core-hole 
site is relevant to the former while the latter vanishes. 

The first term of Eq. (|2.4j) gives the spectral shape as 
a function of uj similar to the conventional one-spin cor- 
relation function familiar to INS theory; the presence of 
/a (w») and /j (w,-) modifies the q c dependence of the 
spectral intensity. The complicated q c dependence is a 
direct consequence of the inclusion of non-local spin ex- 
citation around the core-hole site as seen from Eq. (|2.7|l , 
which is missing in the fast collision approximation. The 
second term of Eq. (|2.4|) gives the spectral shape arising 
from the non-local exchange type excitations occurred 
around the core-hole site. This type of contribution is 
missing in the INS spectra and appears only when the 



theory covers uj dependence beyond the fast collision ap- 
proximation. 

We end this section with the explanation how one- and 
two-triplon excitations are manifested in the one- and 
two-spin correlation functions. For simplicity, we do not 
mention n-triplon excitation with n more than three be- 
cause n — 1 and 2 dominate quantitatively. From the 
definitions Eqs. (|2.7[) and (I2.8[) . one- and two-spin cor- 
relation functions are associated with total spin S = 1 
and spin-conserving excitations, respectively. Since one- 
and two-triplon excitations are total spin S — 1 and 
S = 0, 1, 2, respectively, one-spin correlation function in- 
cludes both one- and two-triplon excitations with S = 1 
while two-spin correlation function includes two-triplon 
excitation alone. 

Note that, in the undoped ladder system, multi-triplon 
contributions with different parity do not mix because 
the system is invariant with respect to reflection about 
the centerline of the ladder^ Therefore, one- and two- 
triplon contributions involved in the one-spin correla- 
tion function (a;,; q c , q ai can be separated. That 
is, one- and two-triplon contributions are found in the 
q a ~ ^ and q a = modes, respectively, in which the spin 
excitations in leg I and 2 are summed up in anti-phase, 
and in phase, respectively. 



III. NUMERICAL RESULTS 

In order to evaluate (Wi)'s and (u;i)'s, we must 
prepare the eigenstates and the corresponding energies 
of -f/mag and -ff m ag> the spin Hamiltonian in the interme- 
diate state. In addition, the incident photon energy u>i 
should be specified, which is chosen as the peak position 
of the absorption coefficient. In the following, we shall 
describe how these preparations are made. 



Eigenstates of i/ ma g and ff^ ag 



We consider a system consisting of 2 x 8 spins of 
S = 1/2 with periodic boundary conditions for the initial 
and final states, as shown in Fig. [T] (b). The exchange 
couplings are chosen here as Jii = 186 meV, J± = 124 
meV, and J cyc = 31 meV j 33 ' 42 " — Representing i/ mag by 
a matrix of 12870 x 12870 dimensions in the subspace of 
the z' component of the total spin S* ot = 0, we diagonal- 
ize the Hamiltonian matrix. We obtain the ground state 
energy as e g /(NJ±) = —0.733. The H' ma& is obtained 
from -ffmag by eliminating the spin at the core-hole site. 
Therefore -ff ma g consists of 15 spins, and may be repre- 
sented by a matrix with 6435 x 6435 dimensions in the 
subspace of S£ ot — ±1/2. 



5 




co/J , co/J 



B. Absorption coefficient 

By substituting the eigenvalues and the eigenstates 
evaluated on finite-size cluster into Eq. (A.l) of ref. [24], 
we obtain the L 2 .3-absorption coefficient Aj(u>i) (J = 3/2 
or j = 1/2). Figure [2] shows the calculated Aj(u>i) 
as a function of the incident photon energy. The ori- 
gin of photon energy is set to be u>i = e corc . The 
e coro stands for the energy required to create a 2p core- 
hole in the multiplet j = 1/2 or 3/2 of the (3d) 10 - 
configuration. For Sri4Cu2404i, we take T/J± = 2.4, 
since a typical value of V ~ 300 meV for Cu Zv3-edge 
and Jj_ = 124 meV. The calculated curve turns out to 
be very close to the Lorentzian shape. The peak position 
is slightly shifted from Wi — e C orc, instead, it is around 
Ui = w? ~ e coro + 1.365Ji for r/J ± = 2.4. 



C. Evaluation of the coefficients 

By using the eigenvalues and eigenfunctions on a two- 
leg ladder of 2 x 8 spins, we calculate the coefficients 
for u>{ = cj°. Table Q] lists the calculated values. For 

fji's, |/a | and I/3 | are rather smaller than |/} | with 
T/J± = 2.4, while they become larger with T/J± in- 
creasing. It implies that the effect of magnetic excita- 
tions on neighboring sites increases with decreasing value 
of T. As regards |/ 2 (2) | ~ 2|/ 3 (2) | » |/| 2) | for 

T/J± = 1.0 ~ 2.4. This suggests that the exchange type 
disturbance does not extend beyond the nearest neighbor 
site in both leg and rung directions. 



FIG. 3: (Color online) Spin correlation functions calculated 
on a 2 x 8 spin ladder as a function of the energy loss uj for 
available q c values. Left panel is one-spin correlation function 
with polarization changing channel and right panel is two- 
spin correlation function with polarization conserving chan- 
nel. Lower and upper panels are for q a = and it, respec- 
tively. The (red) thin lines in panels (a) and (c) represent the 
lower and upper boundaries expected from the two-triplon 
continuum.— The (red) thin line in panel (b) is the dispersion 
curve expected from the one-triplon excitation.— 



D. RIXS spectra 

We calculate the correlation functions from Eqs. 
([23)1 and (pTo) . Figure [3] shows q c , q a , u) and 

y( 2 '(w°;g c ,g a ,cj) numerically calculated with T/J± = 
2.4. It should be reminded that Y^(u)^;q c ,0,u)) and 
Y"W q c ,7r,u)) include the contributions from two- and 
one-triplon excitations, respectively. 

In Fig. [3] (a), the energy profile of Y^ q c , 0, lu) 
is plotted for accessible values of q c . The lower and up- 
per boundaries of the two-triplon continuum, shown in 
Schlappa et al.'s paper j£ are also illustrated. The peak 
position of Y^ (cof; q c ,0,iL>) captures well the dispersive 
behavior of the lower boundary. Quantitative discrep- 
ancy can be ascribed to the finite size effect. Since Eq. 
(|2.5[) gives no finite intensity at the zone center q c = 0, 
our theory, unfortunately, cannot observe the spin gap 
directly. However, a little larger cluster may give the 
better extrapolation of the magnitude of the spin gap at 
the zone center. 

In Fig. [3] (b), we can see that the q c dependence of 
the peak position of Y^ q c , n, w) traces one-triplon 
dispersion relation. 7 Although Schlappa et al. fixed their 
experimental configuration such that the scattering vec- 
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TABLE I: Coefficients (ui)'s and f„ (u>i)'s in units of 1/Jx for different values of T/J± evaluated for 2x8 spin ladder. 



r/Jx 




/ 2 (1) K°) 




/i x) K 9 ) 


2.4 
1.6 
1.0 


(0.034, -0.849) 
(0.063, -1.285) 
(0.019, -2.080) 


(0.052, -0.037) 
(0.084, -0.085) 
(0.110, -0.195) 


(0.029, -0.026) 
(0.045, -0.055) 
(0.056, -0.116) 


(0.035, -0.026) 
(0.056, -0.060) 
(0.073, -0.134) 




r/Jx 










2.4 
1.6 
1.0 




(0.145, -0.108) 
(0.231, -0.246) 
(0.299, -0.550) 


(0.076, -0.060) 
(0.120, -0.113) 
(0.154, -0.292) 


(0.001, 0.004) 
(0.003, 0.004) 
(0.007, 0.003) 



yl'l Y ,2> 




FIG. 4: (Color online) Frequency-integrated intensities 
of the correlation functions (a) I fl \uJi; q c , q a )/N and (b) 
J^ 2 '(cJi ; q c ,q a )/N calculated on a 2 x 8 spin ladder as a func- 
tion of q c . (Black) circles and (red) triangles are intensities 
for q a = and q a = it, respectively. Filled and open symbols 
correspond to r/Jx = 2.4 and 1.0, respectively. Solid and 
broken lines are guides to the eye. 

tor was confined within the fee-plane, RIXS spectra can 
observe the one-triplon dispersion around the zone center 
q c ~ if the projected momentum transfer q a could be 
set as 7T. Such a possibility shall be discussed in Sec. [IV] 

The spectral shape of (u)f; q c , q a , w) as a function 
of ui is shown in Figs. [3] (c) and (d). It seems to have 
more weights at higher uj than those of Y^ l '{u)^\ q c , q a , w). 
This tendency is the same as observed in the analy- 
sis of one-dimensional system^ The peak intensity of 
y( 2 )(w°; q c , 0, lo) is very small when q c — ~ 71-/4, but 
grows rapidly when q c reaches ir/2. Then, the spectrum 
becomes broader and the spectral weight shifts to the 
higher energy region as q c goes to the zone boundary. 
The peak intensity of Y^ 2 \u°; q c , tt, uj) is small at q c = 0, 
peaks at q c — tt/2, and vanishes at the zone boundary. 

Figure @] shows the integrated intensities defined by 

lM(uj?;q e ,q a ) = J V« (u£« ,9 , w)^, (3.1) 
I^(^;q c ,q a ) = jy( 2 >K°;g c , 9a , W )^. (3.2) 




FIG. 5: (Color online) A schematic diagram of the RIXS ex- 
perimental configuration. The incident and scattering angles 
are 9 and tp, respectively. The scattering plane is perpen- 
dicular to the ca-plane. The angle between the scattering 
plane and the c-axis is ip, which is fixed to zero in the actual 
experiment .£ Scattering vector is defined as — q = — q 4 . 



The lW(w$;q e , 0) vanishes with q c — » 0, and increases 
gradually with q c — > tt. The I^(w^;q c ,Tx), on the other 
hand, remains finite at q c = 0, and increases rapidly with 
q c — > tt. The presence of /a (w,) and /j (<*;*) makes 
the g c -dependence deviate from that of the dynamical 
structure factor predicted by the fast collision approx- 
imation. The deviation is small for T/J± = 2.4, but 
becomes conspicuous with T/J± = 1.0, because of the 

increase of and \fP(u>?)\. The /( 2 )(^°; gc ,0) 

is very small but finite at q c = 0, increases as q c grows 
peaking around q c — 37r/4, then remains finite at the 
zone boundary. The l( 2 \u>®; q c , tt) starts finite at q c = 0, 
peaks around q c = tt/2, then vanishes at the zone bound- 
ary. The l( 2 '(u>9; q c ,q a ) is found one order of magnitude 
smaller than I^ 1 ' (u>®; q c , q a ) over the entire Brillouin zone. 



IV. APPLICATION TO Sri 4 Cu240 4 i 

Now we attempt to compare our results with those 
observed by the RIXS experiment in two-leg spin lad- 
der system Sri4Cu2404if£ The legs and rungs of the lad- 
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FIG. 6: (Color online) RIXS spectra as a function of the en- 
ergy loss uj for the available momentum transfer projected 
on the c-axis (q c ) in the a polarization. Solid (black) and 
long-dotted (red) lines show the total intensity and the contri- 
butions from Y ' 2 ' (ojf; q c , q a , ui), respectively. The calculated 
curves are convoluted with the Lorentz function with the half- 
width of half-maximum 78 meV. (a) For q a — 0. Filled circles 
are the experimental datai Thin (blue) curves are the lower 
and upper boundaries expected from two-triplon continuum. 7 
(b) For — q a = tt. Thin (blue) line is the dispersion curve ex- 
pected from one-triplon excitation.— 



ders are along the c and a axes with lattice constants 
a = 11.459A, b = 13.368A, and along the c-direction the 
unit cell for the ladders is cl = 3.93lAj22 The lattice 
parameter is a r — 3.84A along the a-directions for the 
ladders M 



A. L3-edge spectrum 

In this section, we compare our results with the ex- 
periment observed at the Cu L3-edge in Sri4Cu2404i. 
Figure [5] illustrates a schematic sketch of the experimen- 
tal geometry performed by Schlappa et al. 7 Notice that 
the direction of the scattering vector in the experiment is 
opposite to that of ours. In the following, when the finite 



numerical values of the momentum transfer q c and q a are 
mentioned, we adopt the experimental definition. The 
surface of the sample is perpendicular to the 6-axis. Since 
the scattering plane is fixed parallel to the bc(zx) -plane, 
the momentum transfer along the rung direction is zero 
(<Za = 0). Hence, the experiment probed Y^'(u)i; q c , 0,uj) 
and Y^ 2 '(uii;q c ,0,uj). The incident and scattering an- 
gles are 9 and ip, respectively. In the experiment, both 
ip = 90° and 130° were used. We choose the latter since 
it covers approximately 90 % of the Brillouin zone along 
the c-direction at the Cu L3-edge. By changing 9, the 
RIXS spectra for different q c were obtained. The polar- 
ization vector of the incident photon is then expressed 
in the xyz coordinate as oti = (0,-1,0) for the a po- 
larization and on — (xfjOiXD f° r t ne 71 polarization. 
Similarly, the polarization of the scattered photon is ex- 
pressed as af = (0,-1,0) for the a' polarization and 
atf = (x/jOiX/) f° r the it' polarization. The polariza- 
tion is usually separated with the incident photon, but 
not separated with the scattered photon in experiments. 
In such a situation, we may express the RIXS spectra 
depending on the polarization of the incident photon as 

I(ui;q c ,q a ,uj) 



(^pj Y^\ujv,q c , qa, u) + Y^(u i; q c , q a ,w) 



(4.1) 



for cr-polarization, and 

I(u!i;q c ,q a ,uj) 



X, 



YW (cu t ; q c , q a , u) + {x}xl?Y {2) q c , q a , u) 



(4.2) 



for 7r-polarization. The contribution of Y^ 2 \uii, q c ,q a , uj) 
relative to that of Y^ (w,i; q c , q a , ui) is enhanced by 
(2/xp 2 in the a polarization. The contribution of 
Y( 2 ) (uJi, q c , q a , uj) in the tt polarization is smaller than 
that in the a polarization by a factor (xjxf) 2 - 

Figures|5](a) and (b) show the RIXS spectra as a func- 
tion of energy loss uj with — q a = and 7r, respectively, in 
the cr polarization. We put T/J± = 2.4. The calculated 
curves are convoluted by the Lorentzian function with 
the half-width-half-maximum of the possible resolution, 
78 meV. Since the contribution from Y^luji, q c , q a , uj) is 
larger than that from Y^ {uji, q c , q a , uj) in most condi- 
tions, the dispersive behavior of the lower boundary of 
two-triplon continuum as well as one-triplon dispersion 
can be well reproduced as is shown in Fig. [3] 

Since the experiment has measured q a — mode alone, 
the one-triplon dispersion explained above has not de- 
tected yet. Here, let us discuss a possible configuration 
to access —q a = n mode by rotating sample by angle tp 
around the 6-axis (see Fig. [5]). An easy way to access 
— q a = 7r mode is setting cp = 90°. Then, we can see that 
when ip > 138° at the Cu L3-edge, the incident angle 9 
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FIG. 7: (Color online) RIXS spectra at the Cu M3-edge as 
a function of the energy loss uj for the available momentum 
transfer projected on the c-axis (q c ) in the a polarization. 
Solid (black) and long-dotted (red) lines show the total in- 
tensity and the contributions from (cj°; q c , q a , uj), respec- 
tively. The calculated curves are convoluted with the Lorentz 
function with the half-width of half-maximum 78 meV. 



becomes positive, i.e., such a configuration is mathemat- 
ically allowed. Examples are 9 ~ 2.5° for ip = 140°, and 
9 ~ 9° for ip = 150°. As seen from Figs. M (a) and (b), 
we find that the intensity expected from —q a = ^ mode 
is roughly the same order of magnitude observed in the 
experiment for q a = 0. Therefore, we assert that local 
minimum of the one-triplon dispersion at — q =—(q cl Q a ) 
=(0,7r) can be observed by the RIXS spectrum. The 
larger ip can be, the larger attainable maximum value of 
q c becomes away from the zone center. 

Next, we concentrate on the contribution of 
y( 2 )(w°; q c , q a , uj), which is smaller than that of 
(Uj ; q c ,q a ,u>) even in the presence of the geometrical 
coefficients appeared in Eqs. (|4.1[) and (|4.2p . Though it 
is small, the contribution of Y^> (w?; q c , q a , uj) has weight 
in the higher energy transfer region. The intensity grows 
larger as \q c \ increases with \q c \ — > 7r/2, then decreases 
toward the zone boundary. In Schlappa et aVs data, for 
each q c , the intensity is accumulated within the region 
centered at the peak position with half- width about 0.1 
eV. However, around — q c = — 0.8n ~ — 0.6-7T, the tail of 
the intensity extends to the higher energy side about 0.2 
eV. This trend coincides with that shown by the two-spin 
correlation function, which leads us to speculate that the 
observed weak intensities in the high energy region may 
be attributed to the contribution of Y ^ (cj°; q c ,q a: uj). 



B. A/3-edge spectrum 

Our analysis has shown that the contribution of 
Y( 2 \u)®; q C: q a , uj) plays a minor role in Sri4Cu2404i in 



making a quantitative description of the L-edge pro- 
cesses. On the other hand, we have confirmed that 
when the ratio Y / J±_ becomes smaller, the contribution 
of Y( 2 \ui®] q c , q a , uj) becomes larger. Therefore, we as- 
sert that the RIXS signal at the Cu M2 i 3-edge may con- 
tain significant contribution of Y^ (w?; q c , q a: oj). Let us 
study such a situation and examine the possibility. 

The present theory is applicable to analyze the mag- 
netic excitation spectrum at the Cu M2,3-edge. In the 
following, we focus on the M3-edge case. A prominent 
difference between the Ma-edge and Z^-edge events ap- 
pears in the values of parameters. Since the core-hole 
is in the 3p instead of 2p state, the resonant energy at 
the M3-edge (747 eV) becomes smaller than that at the 
L3-edge (931 eV). Thus, the core-hole lime-time broad- 
ening width becomes smaller value, for instance, T = 0.2 
eV^£ which results in T/J± = 1.6. The smaller reso- 
nant energy leads us to the different value of IqJ. This 
makes geometrical condition the equipment should sat- 
isfy tougher, but attainable. Figure [7] shows the calcu- 
lated result of the M3-edge spectrum for q a = with 
ip = 130° and a incident polarization. In this case, 
the transferred momentum projected along the leg (q c ) 
can be covered up to nearly 80% of the first Brillouin 
zone. The contribution of the two-spin correlation func- 
tion to the total spectrum at the M3-edge is larger than 
that at the L3-edge. However, it is hard to distin- 
guish pure Y^(uj^;q c ,0,uj) contribution since it is hid- 
den by the large Y^iuif; q c , 0, uj) contribution for most 
q c . An important difference is found at — q c = 0.757T 
where the total intensity is almost exclusively generated 
by Y^(uj^;q c ,0,uj). This value of q c is near the geo- 
metrical limit but within the reach. Experiments have 
been reported in other Cu-oxides such as CaCu02 and 
SrCu0 2 at the Cu Af 2 , 3 -edges4^ We hope the RIXS 
measurement of at Cu M 2 ,3-edge in Sri4Cu2404i will be 
carried out in near future. 



V. SUMMARY AND DISCUSSION 

We have studied the magnetic excitation spectra of the 
L-edge RIXS in undoped cuprate, in particular, in quasi 
one-dimensional two-leg spin ladder system. We have an- 
alyzed the second-order dipole allowed process through 
the intermediate state, in which there is no spin degree 
of freedom at the core-hole site. This nature of the in- 
termediate state is found to affect strongly the transition 
amplitudes of spin excitations. Then, the RIXS spectra 
have been derived as the one-spin and two-spin correla- 
tion functions in the channels with and without chang- 
ing polarization, respectively. Note that they include the 
contributions of the magnetic excitations not only on the 
core-hole site but also on the neighboring sites. The cor- 
relation functions are expressed by a set of numerical 
coefficients reflecting the weight of each magnetic excita- 
tions. 

Once the expressions of the correlation functions are 
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obtained, we could concentrate on evaluating them with 
approximation methods available. The coefficients could 
be rather accurately obtained numerically in a system 
with small size, since the relevant excitations are re- 
stricted around the core-hole site. We have evaluated 
them and, subsequently, the correlation functions in a 
finite-size two-leg ladder. In the two-leg ladder system, 
the one-spin correlation function includes both one- and 
two-triplon excitations with S = 1. We find one- and 
two-triplon excitations emerge separately by choosing 
rung wave number — q a = ix and 0, respectively. An 
application to Sri4Cu2404i has revealed that the cal- 
culated RIXS spectrum captures well the dispersive be- 
havior shown by the lower boundary of two-triplon con- 
tinuum. By adjusting the geometrical configuration of 
the experiment, one-triplon dispersion around the zone 
center could be detectable in the RIXS measurement. 
The observed weak intensity in the higher energy region 
around — q c = — 0.8n ~ — 0.67T might be the contribution 
from the two-spin correlation function, which could be 
clearly detected at the Cu M 3 -edge. 

We finally comment on the plausibility to adopt spin 
only model to explain the current experiment. Although 
our theory has reproduced semi-quantitatively well the 
lower boundary of the two-triplon continuum over a wide 



range of the Brillouin zone, the corresponding spin cor- 
relation function has no intensity at the zone center. On 
the other hand, a small cluster analysis on the Hubbard 
model presented by Schlappa et al. has given a finite in- 
tensity at q c = 0. Since the RIXS processes in the actual 
materials are very complicated, there is a possibility that 
the observed RIXS signals involves something missing in 
the spin correlation function. For the if-edge RIXS spec- 
tra, for instance, Jia et al. has demonstrated that the 
RIXS intensity and the spin dynamical structure factor 
show difference in several systems 4£ Whether the same 
is true to the present system is an intriguing problem and 
relegated to a future study. 
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